LIE GROUPS AND MECHANICS, 
AN INTRODUCTION 



BORIS KOLEV 



Abstract. The aim of this paper is to present aspects of the use of Lie 
groups in mechanics. We start with the motion of the rigid body for 
which the main concepts are extracted. In a second part, we extend the 
theory for an arbitrary Lie group and in a third section we apply these 
methods for the diffeomorphism group of the circle with two particular 
examples: the Burger equation and the Camassa-Holm equation. 



Introduction 

The aim of this article is to present aspects of the use of Lie groups in 
mechanics. In a famous article Arnold showed that the motion of the 
rigid body and the motion of an incompressible, inviscid fluid have the same 
structure. Both correspond to the geodesic flow of a one-sided invariant 
metric on a Lie group. From a rather different point of view, Jean-Marie 
Souriau has pointed out in the seventies [22] the fundamental role played 
by Lie groups in mechanics and especially by the dual space of the Lie 
algebra of the group and the coadjoint action. We aim to discuss some 
aspects of these notions through examples in finite and infinite dimension. 
The article is divided in three parts. In Section ^ we study in detail the 
motion of an n-dimensional rigid body. In the second section, we treat the 
geodesic flow of left-invariant metrics on an arbitrary Lie group (of finite 
dimension). This permits us to extract the abstract structure from the case 
of the motion of the rigid body which we presented in Section ^ Finally, 
in the last section, we study the geodesic flow of H'' right-invariant metrics 
on Diff(S^), the diffeomorphism group of the circle, using the approach 
developed in Section [2 Two values of k have significant physical meaning 
in this example: A: = corresponds to the inviscid Burgers equation jl6j and 
k = 1 corresponds to the Camassa-Holm equation 01^. 



1. The motion of the rigid body 

1.1. Rigid body. In classical mechanics, a material system (E) in the am- 
bient space is described by a positive measure fi on with compact 
support. This measure is called the mass distribution of (S). 

• If is proportional to the Dirac measure 6p, (S) is the massive point 
P, the multiplicative factor being the mass m of the point. 

• If /X is absolutely continuous with respect to the Lebesgue measure 
A on M'^, then the Radon-Nikodym derivative of ^ with respect to A 
is the mass density of the system (S). 
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In the Lagrangian formalism of Mechanics, a motion of a material system 
is described by a smooth path of embeddings of the reference state T, = 
Supp(fj,) in the ambient space. A material system (S) is rigid if each map 
If is the restriction to S of an isometry g of the Euclidean space M^. Such a 
condition defines what one calls a constitutive law of motion which restricts 
the space of probable motions to that of admissible ones. 

In the following section, we are going to study the motions of a rigid 
body (S) such that S = Supp{fi) spans the 3 space. In that case, the 
manifold of all possible configurations of (S) is completely described by the 
6-dimensional bundle of frames of M^, which we denote 7^(M^). The group 
of orientation-preserving isometrics of acts simply and transitively 
on that space and we can identify 7^(M^) with D3. Notice, however, that 
this identification is not canonical - it depends of the choice of a "reference" 
frame Kq. 

Although the physically meaningful rigid body mechanics is in dimension 
3, we will not use this peculiarity in order to distinguish easier the main 
underlying concepts. Hence, in what follows, we will study the motion of an 
re-dimensional rigid body. 

Moreover, since we want to insist on concepts rather than struggle with 
heavy computations, we will restrain our study to motions of a rigid body 
having a fixed point. This reduction can be justified physically by the possi- 
bility to describe the motion of an isolated body in an inertial frame around 
its center of mass. In these circumstances, the configuration space reduces 
to the group SO{n) of isometrics which fix a point. 

1.2. Lie algebra of the rotation group. The Lie algebra 5o(re) of SO{n) 
is the space of all skew-symmetric n x re matrices^. There is a canonical 
inner product, the so-called Killing form |2S] 

which permit us to identify so (re) with its dual space 50 (n)*. 
For X and y in M", we define 

L*(x, y){9) = (Qx) -y, 17 G 5o(re) 

which is skew-symmetric in x,y and defines thus a linear map 

2 

L* : /\M" ^so(re)* . 

This map is injective and is therefore an isomorphism between so (re)* and 
/\ M", which have the same dimension. Using the identification of so(n)* 
with so(n), we check that the element L{x,y) of so(n) corresponding to 
L*{x,y) is the matrix 

(1) L{x,y) = yx^ - xy^ . 

where x* stands for the transpose of the column vector x. 



^In dimension 3, we generally identify the Lie algebra so (3) with R'^ endowed with the 
Lie bracket given by the cross product lui x uj2- 
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1.3. Kinematics. The location of a point a of the body S is described by 
the column vector r of its coordinates in the frame Kq- At time t, this point 
occupies a new position r{t) in space and we have r{t) = g{t)r, where g{t) is 
an element of the group S0{3). In the Lagrangian formalism, the velocity 
v(a, t) of point a of S at time t is given by 

d 

v(a,t) = -Q^^ia,t) = g{t) r. 
The kinetic energy K of the body S at time t is defined by 

(2) K(t) = i^||v(a,t)f df, = ^ljgrf df, = ^ JjQrf df, 

where Q, = g^^ g lies in the Lie algebra so(n). 

Lemma 1.1. We have K = — ^ ti{Vljn), where J is the symmetric matrix 
with entries 

r 

Jy. 

Proof. Let L : /\^M"' so(n) be the operator defined by (0). We have 

(3) L{r,9.r) = {rr*)n + 9.{rr*) 1^ G 5o(n), r G M", 
where rr* is the symmetric matrix with entries XiXj. Therefore 

(4) {Vlr) ■{Vlr)= L*{r,nr)n = tr (L(r,l^r)l7) = - tr (l7(rr*)17) , 

which leads to the claimed result after integration. □ 

The kinetic energy K is therefore a positive quadratic form on the Lie 
algebra so (n). A linear operator A : so(n) 5o{n), called the inertia tensor 
or the inertia operator, is associated to K by means of the relation 

K = ^{A{n),n) , JlGso(n). 

More precisely, this operator is given by 

(5) A{n) = jn + nj = [ {n rr^ + rr^n) dfi . 

Remark. In dimension 3 the identification between a skew-symmetric matrix 
Q and a vector oj is given by ivi = —^223) = ^13 and UJ3 = —i^i2- If we 
look for a symmetric matrix / such that A{Q) correspond to the vector lu, 
we find that 

+ —xy —xz \ 

—xy + z"^ —yz dfi, 

—xz —yz x^ + y"^ J 

which gives the formula used in Classical Mechanics. 
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1.4. Angular momentum. In classical mechanics, we define the angular 
momentum of the body as the following 2-vector^ 

M{t)= / (gr) A (gr) d^i . 

Lemma 1.2. We have L{M) = gA{Q)g^^. 
Proof. A straightforward computation shows that 

L{gr,gr) = gQrr^g~^ + grr''i}g~^ . 

Hence 

L{M) = J L{gr,gr) d^i = gA{n)g~^ . 

□ 

1.5. Equation of motion. If there are no external actions on the body, 
the spatial angular momentum is a constant of the motion. 

Coupled with the relation L{M.) = gA{Q)g ^, we deduce that 

(7) A{tl) = A{^)Vl - VlA{Vt) 

which is the generalization in n dimensions of the traditional Euler equation. 
Notice that if we let M = A{Q), this equation can be rewritten as 

(8) M=[M,rL]. 

1.6. Integrability. Equation has the peculiarity that the eigenvalues 
of the matrix M are preserved in time. Usually, integrals of motion help to 
integrate a differential equation. The Lax pairs technique jEl is a method 
to generate such integrals. Let us summarize briefly this technique for finite 
dimensional vector spaces. Let ii = F{u) be an ordinary differential equation 
in a vector space E. Suppose that we were able to find a smooth map 
L : E ^ End(F), where F is another vector space of finite dimension, 
with the following property: if u{t) is a solution of tt = F{u), then the 
operators L{t) = L{u{t)) remain conjugate with each other, that is, there is 
a one-parameter family of invertible operators P{t) such that 

(9) m = p{tr'mp{t) . 

In that case, differentiating Q, we get 

(10) L=[L,B] 

where B = P^^P. Conversely, if we can find a smooth one-parameter family 
of matrices B{t) € End{F), solutions of equation ((TIH) . then © is satisfied 
with P{t) a solution of P = PB . If this is the case, then the eigenvalues, 
the trace and more generally all conjugacy invariants of L{u) constitute a 
set of integrals for ii = F{u). 
A Hamiltonian system on R2JV 

called completely integrable if it has N 
integrals in involution that are functionally independent almost everywhere. 

^In the Euclidean 3-space, 2-vectors and 1-vectors coincide. This is why, usually, one 
consider the angular momentum as a 1-vector. 
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A theorem of Liouville describes in that case, at least quahtatively, the 
dynamics of the equation. This is the reason why it is so important to find 
integrals of motions of a given differential equation. 

Using the Lax pairs technique, Manakov [2^ proved the following theorem 

Theorem 1.3. Given any n, equation (jH)) has 



N{n) = i 



n(n 



integrals of motion in involution. The equation of motion of an n-dimensional 
rigid body is completely integrable. 

Sketch of proof. The proof is based on the following basic lemma. 

Lemma 1.4. Euler's equations (jH)) of the dynamics of an n-dimensional 
rigid body have, for any n, a representation in Lax's form in matrices, 
linearly dependent on a parameter A € C, given by L\ = M + J^X and 

Bx = n + jx. 

Hence, the polynomials = tr [M + J'^X)^, {k = 2, . . . ,n) are time- 

independent and the coefficients Pk{X) are integrals of motion. Since M is 
skew-symmetric and J is symmetric, the coefficient of A* in Pfc(A) is nonzero, 
provided s has the same parity as k. The calculation of N{n) here presents 
no difficulties. □ 



2. Geodesic flow on a Lie Group 

In this section, we are going to study the geodesic flow of a left invariant 
metric on a Lie group of finite dimension. Our aim is to show that all the 
computations performed in Section ^ are a very special case of the theory 
of one-sided invariant metrics on a Lie group. Later on, we will use these 
techniques to handle partial differential equations. We refer to [2] from 
where materials of this section come from and to Souriau's book [25! 
a thorough discussion of the role played by the dual of the Lie algebra in 
mechanics and physics. 

2.1. Lie Groups. A Lie group G is a group together with a smooth struc- 
ture such that g i— > g^^ and {g, h) i— >■ gh are smooth. On G, we define 
the right translations Rh : G ^ G hy Rg{h) = hg and the left translations 
Lg-.G^Ghy Lg{h) = gh. 

A Lie group is equipped with a canonical vector-valued one form, the so 
called Maurer-Cartan form Lo{Xg) = Lg~iXg which shows that the tangent 
bundle to G is trivial TG ~ G x g. Here g is the tangent space at the group 
unity e. 

A left- invariant tensor is completely defined by its value at the group 
unity e. In particular, there is an isomorphism between the tangent space 
at the origin and left-invariant vector fields. Since the Lie bracket of such 
fields is again a left-invariant vector field, the Lie algebra structure on vector 
fields is inherited by the tangent space at the origin g. This space g is called 
the Lie algebra of the group G. 
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Remark. One could have defined the Lie bracket on g by pulhng back the 
Lie bracket of vector fields by right translation. The two definitions differ 
just by a minus sign 

Example. The Lie algebra so(n) of the rotation group SO{n) consists of 
skew-symmetric n x n matrices. 

2.2. Adjoint representation of G. The composition Ig = Rg-iLg : G — > 
G which sends any group element h E G to ghg~^ is an automorphism, that 
is, 

Ig{hk) = Ig{h)lg{k). 

It is called an inner automorphism of G. Notice that Ig preserves the group 
unity. 

The differential of the inner automorphism Ig at the group unity e is 
called the group adjoint operator Adg defined by 

Adg-.g^Q, Adg u = ■^\t=o Ig{h{t)), 

where h{t) is a curve on the group G such that /i(0) = e and /i(0) = lo G g = 
Tf,G. The orbit of a point a; of g under the action of the adjoint representa- 
tion is called an adjoint orbit. The adjoint operators form a representation 
of the group G (i.e. Adgh = AdgAdh) which preserves the Lie bracket of q, 
that is, 

[Adg^,AdgLO] =Adg , LO ] . 

This is the Adjoint representation of G into its Lie algebra 0. 
Example. For g G SO{n) and € 5o(n), we have Adg = gQ,g~^ . 

2.3. Adjoint representation of q. The map Ad, which associates the 
operator Adg to a group element g £ G, may be regarded as a map from 
the group G to the space End{g) of endomorphisms of g. The differential 
of the map Ad at the group unity is called the adjoint representation of the 
Lie algebra g into itself, 

ad: g^ End{g), ad^ u; = ■^\t=o Adg(^t) u). 

Here g{t) is a curve on the group G such that g{0) = e and g{0) = ^. Notice 
that the space {ad^ u, e g} is the tangent space to the adjoint orbit of the 
point Lo E g. 

Example. On the rotation group SO{n), we have ads ^ = [S;^^], where 
[H, ] = 30, — is the commutator of the skew-symmetric matrices S and 
Q. As we already noticed, for n = 3, the vector [(,,uj] is the ordinary cross 
product ^ X w of the angular velocity vectors ^ and a; in M^. More generally, 
if G is an arbitrary Lie group and [!^,uj] is the Lie bracket on g defined 
earlier, we have ad^ a; = [^,a;]. 
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2.4. Coadjoint representation of G. Let g* be the dual vector space to 
the Lie algebra g. Elements of g* are linear functionals on g. As we shall see, 
the leading part in mechanics is not played by the Lie algebra itself but by 
its dual space g*. Souriau [251 pointed out the importance of this space in 
physics and called the elements of g* torsors of the group G. This definition 
is justified by the fact that torsors of the usual group of affine Euclidean 
isometrics of represent the torsors or torques of mechanicians. 

Let A : E ^ F he a. linear mapping between vector spaces. The dual 
(or adjoint) operator A* , acting in the reverse direction between the corre- 
sponding dual spaces. A* : F* ^ E* , is defined by 

{A* a){x) = a{A x) 

for every x £ E, a F*. 

The coadjoint representation of a Lie group G in the space g* is the repre- 
sentation that associates to each group element g the linear transformation 

Ad; : g* ^ g* 

given by Ad; = {Adg-i)*. In other words, 

{Ad; m){io) = m{Adg u) 

for every g £ G, m £ q* and a; € g. The choice of g^^ in the definition of 
Ad; is to ensure that Ad* is a left representation, that is Ad;/^ = Ad;Ad^ 
and not the converse (or right representation). The orhit of a point m of g* 
under the action of the coadjoint representation is called a coadjoint orhit. 
The Killing form on g is defined by 

A;(,^, uj) = tr (adg adiS) ■ 

Notice that k is invariant under the adjoint representation of G. The Lie 
group G is semi-simple if k is non-degenerate. In that case, k induces an 
isomorphism between g and g* which permutes the adjoint and coadjoint 
representation. The adjoint and coadjoint representation of a semi-simple 
Lie group are equivalent. 

Example. For the group 50(3) the coadjoint orbits are the sphere centered 
at the origin of the 3-dimensional space so (3)*. They are similar to the 
adjoint orbits of this group, which are spheres in the space so (3). <0 

Example. For the group SO{n) (n > 3), the adjoint representation and 
coadjoint representations are equivalent due to the non-degeneracy of the 
Killing form^ 

A:(H,J7) = itr (Hl^*), 

where il* is the transpose of O relative to the corresponding inner product 
ofM". Therefore 

Ad;M = gMg-\ 
for M £ 50{n)* and g € SO{n). 



This formula is exact up to a scaling factor since a precise computation for so{n) gives 
k{X,Y) = (7i-2)tr(XY). 
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Despite the previous two examples, in general the coadjoint and the ad- 
joint representations are not alike. For example, this is the case for the 
Poincare group (the non-homogenous Lorentz group) cf. |13j . 

2.5. Coadjoint representation of q. Similar to the adjoint representation 
of g, there is the coadjoint representation of q. This later is defined as the 
dual of the adjoint representation of g, that is, 

ad* : g —i- End{Q*), ad^ m = {ad^)*{m) = — — \f^Q ^'ig(t) 

where g{t) is a curve on the group G such that ^(0) = e and g{0) = S^. 

Example. For U e so(n) and M e 5o{n)*, we have ad^ M = — M]. 

Given m € g*, the vectors m, with various C G S; constitute the 
tangent space to the coadjoint orbit of the point m. 

2.6. Left invariant metric on G. A Riemannian or pseudo-Riemannian 
metric on a Lie group G is left invariant if it is preserved under every left 
shift Lg, that is, 

{Xg,Yg)^ = {L,Xg,Lf,Yg),^, g,hGG. 

A left-invariant metric is uniquely defined by its restriction to the tangent 
space to the group at the unity, hence by a quadratic form on g. To such a 
quadratic form on g, a symmetric operator ^ : g — > g* defined by 

{^,uj) = {AC,oj) = {Auj,^), ^,ojeQ, 

is naturally associated, and conversely^. The operator A is called the inertia 
operator. A can be extended to a left-invariant tensor Ag : TgG TgG* 
defined by Ag = L*_iALg~i. More precisely, we have 

{X, Y)g = {AgX, Y)^ = [AgY, X )^ , X, y G TgG. 

The Levi-Civita connection of a left-invariant metric is itself left-invariant: 
if La and Lf, are left-invariant vector fields, so is V l^L^,. We can write down 
an expression for this connection using the operator B : g x g — s- g defined 

by 

(11) {[a,h],c) = {B{c,a),h) 
for every a,b,c in g. An exact expression for B is 

B{a,b) = A-^adl{Aa). 
With these definitions, we get 

(12) iVLM){e) = \[a,b]- ^{Bia, b) + B{b, a)} 



The round brackets correspond to the natural pairing between elements of g and g*. 



LIE GROUPS AND MECHANICS 9 

2.7. Geodesies. Geodesies are defined as extremals of the Lagrangian 
(13) £{g) = J K{g{t),g{t)) dt 

where 

is called the kinetic energy or energy functional. 

If g{t) is a geodesic, the velocity g{t) can be translated to the identity via 
left or right shifts and we obtain two elements of the Lie algebra g, 

called the left angular velocity, respectively the right angular velocity. Let- 
ting m = Ag g & TgG*, we define the left angular momentum tul and the 
right angular momentum by 

niL = L*gm G g*, tur = RgW, G Q*. 

Between these four elements, we have the relations 

(15) LOR = AdgUJL, mR = Ad*gmL, m,L = AujL- 
Note that the kinetic energy is given by the formula 

(16) K = ]^{g,g)g = ]^ {lvl, ^l) = ^ {rriL, ^l) = ^ {Ag g,g)g . 
Example. The kinetic energy of an n-dimensional rigid body, defined by 

(17) K{t) = ^ljgrf di^ = -^tv{njn) 

is clearly a left-invariant Riemannian metric on SO{n). In this example, we 
have = i^L and M = mi. Physically, the left-invariance is justified by 
the fact that the physics of the problem must not depend on a particular 
choice of reference frame used to describe it. It is a special case of Galilean 
invar iance. <) 

2.8. Euler- Arnold equation. The invariance of the energy with respect 
to left translations leads to the existence of a momentum map fi : TG — > q* 
defined by 

dK 

Ki9,m0 = ^^i = (9,RgOg = {m,Rgi) = {R;m,^) = mR(0, 

where is the right-invariant vector field generated by ^ G g. According 
to Noether's theorem |25) . this map is constant along a geodesic, that is 



(18) ^ = 0. 
^ ^ dt 

As we did in the special case of the group SO{n), using the relation 
Ad*g tul and computing the time derivative, we obtain 

(19) —^=ad^^mL. 
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(21) 



This equation is known as the Arnold-Euler equation. Using uji = A ^ rriL, 
it can be rewritten as an evolution equation on the Lie algebra 

(20) ^ = B{u;l,u;l). 

Remark. The Euler- Lagrange equations of problem (jlZj) are given by 

g = LgUOL, 

LOL =B{uJl,UJl)- 

If the metric is bi-invariant, then B{a,b) = for all a, 6 € and ujl is 
constant. In that special case, geodesies are one-parameter subgroups, as 
expected. 

2.9. Lie-Poisson structure on q*. A Poisson structure on a manifold M 
is a skew-symmetric bilinear function {, } that associates to a pair of smooth 
functions on the manifold a third function, and which satisfies the Jacobi 
identity 

{{f,g},h} + {{g,h},f} + {{h,f},g} = 

as well as the Leibniz identity 

{f,9h} = {f,g}h + g{f,h}. 

On the torsor space q* of a Lie group G, there is a natural Poisson struc- 
ture defined by 

(22) {f,9}{m) = {m,[dmf,dmg]) 

for m € g* and f,g& C°°{g*). Note that the differential of / at each point 
m G g* is an element of the Lie algebra g itself. Hence, the commutator 
[dm f,dmg] is also a vector of this Lie algebra. The operation defined above 
is called the natural Lie-Poisson structure on the dual space to a Lie algebra. 
For more materials on Poisson structures, we refer to [20n26j . 

Remark. A Poisson structure on a vector space E is linear if the Poisson 
bracket of two linear functions is itself a linear function. This property is 
satisfied by the Lie-Poisson bracket on the torsors space g* of a Lie group 
G. 

To each function H on a Poisson manifold M one can associate a vector 
field defined by 

L^^f = {HJ} 
and called the Hamiltonian field of H. Notice that 

[^F,Ch] = C{F,H}- 

Conversely, a vector field v on a Poisson manifold is said to be Hamiltonian 
if there exists a function H such that v = £,h- 

Example. On the torsors space g* of a Lie group G, the Hamiltonian field 
of a function H for the natural Lie-Poisson structure is given by S,h{^) = 

Let A be the inertia operator associated to a left-invariant metric on G. 
Then equation (|19|) on g* is Hamiltonian with quadratic Hamiltonian 

H{m,) = ^ [A^^ m,M) , meg*. 
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which is nothing else but the kinetic energy expressed in terms of m = Alo. 
Notice that since ?n.i(t) = Ad*(^^^ nT-R where mji G g* is a constant, each 
integral curve mL{t) of this equation stays on a coadjoint orbit. 

A Poisson structure on a manifold M is non-degenerate if it derives from 
a symplectic structure on M. That is 

{f,9} = , 

where w is a non-degenerate closed two form on M. Unfortunately, the Lie- 
Poisson structure on g* is degenerate in general. However, the restriction 
of this structure on each coadjoint orbit is non-degenerate. The symplectic 
structure on each coadjoint orbit is known as the Kirillov^ form. It is given 

by 

oj{ad*^ m, adl m) = (m, [a,b]) 

where a,6 € g and m € g*. Recall that the tangent space to the coadjoint 
orbit of m G g* is spanned by the vectors a(i| m where ^ describes g. 

3. Right-invariant metric on the diffeomorphism group 

In ^ , Arnold showed that Euler equations of an incompressible fluid may 
be viewed as the geodesic flow of a right-invariant metric on the group of 
volume-preserving diffeomorphism of a 3-dimensional Riemannian manifold 
M (filled by the fluid). More precisely, let G = Dif f^{M) be the group 
of diffeomorphisms preserving a volume form on some closed Riemannian 
manifold M. According to the Action Principle, motions of an ideal (incom- 
pressible and inviscid) fluid in M are geodesies of a right-invariant metric 
on Dif f^{M). Such a metric is defined by a quadratic form K (the kinetic 
energy) on the Lie algebra X^{M) of divergence- free vector fields 

K = \ I M^d^i 

^ JM 

M 2 

where \\v\\ is the square of the Riemannian length of a vector field v G 
X{M). An operator B on ^^(M) x ^^(M) defined by the relation 

{[u,v] = {B{w, u),v) 

exists. It is given by the formula 

B{u, v) = curlii X t; -|- gradp , 

where x is the cross product and p a function on M defined uniquely (mod- 
ulo an additive constant) by the condition divi? = and the tangency of 
B(u, v) to dM. The Euler equation for ideal hydrodynamics is the evolution 
equation 

du 

(23) — = u X curlti — gradp. 

If at least formally, the theory works as well in infinite dimension and the 
unifying concepts it brings form a beautiful piece of mathematics, the details 
of the theory are far from being as clear as in finite dimension. The main 
reason of these difficulties is the fact that the diffeomorphism group is just a 



Jean-Marie Souriau has generalized this construction for other natural G-actions on 
g* when the group G has non null symplectic cohomology |25|. 
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Frechet Lie group, where the main theorems of differential geometry hke the 
Cauchy-Lipschitz theorem and the Inverse function theorem are no longer 
valid. 

In this section, we are going to apply the results of Section to study the 
geodesic flow of a right-invariant metrics on the diffeomorphism group 
of the circle This may appear to be less ambitious than to study the 3- 
dimensional diffeomorphism group. However, we will be able to understand 
in that example some phenomena which may lead to understand why the 
3-dimensional ideal hydrodynamics is so difficult to handle. Moreover, we 
shall give an example where things happen to work well, the Camassa-Holm 
equation. 

3.1. The diffeomorphism group of the circle. The group Diff(S^) is 
an open subset of C°°(§^,§-'^) which is itself a closed subset of C°°(S^,C). 
We define a local chart (C/qj^o) around a point (/?o € Diff{E^) by the 
neighborhood 

Uo = {llv? - V>o\\c"{s^) < 1/2} 

of (fQ and the map 

1 



ZTTl 

The structure described above endows Z)i//(S^) with a smooth manifold 
structure based on the Frechet space C°°(S^). The composition and the in- 
verse are both smooth maps -Di//(S^) x Dif f(S^) — > Diff{S^), respectively 
DiffiS^) Diff{8^), so that L>i//(§^) is a Lie group. 

A tangent vector V at a point (p G Diff{S ) is a function F : ^ TS^ 
such that Tr{V{x)) = ^{x). It is represented by a pair f) € Diff{S^) x 
C°°(§^). Left and right translations are smooth maps and their derivatives 
at a point tp G Diff{S^) are given by 

R^V = {^{i;),vm 
The adjoint action on g = Vect{§^) = C°^(S^) is 

whereas the Lie bracket on the Lie algebra Tjii Dif f{S^) = Vect{S^) = 
C°°(Si) of Diff{S^) is given by 

[u, v] = —{UxV — UVx), U, V € C°°{S^) 

Each V S Vect{S^) gives rise to a one-parameter subgroup of diffeomor- 
phisms •)} obtained by solving 

(24) Vt = v{v) in C^iS') 

with initial data 7]{0) = Id £ Diff{S^). Conversely, each one-parameter 
subgroup t 1-^ 7]{t) G Diff{S^) is determined by its infinitesimal generator 

Evaluating the flow t 1-^ r](t, •) of H24|) at t = 1 we obtain an element expi^{v) 
of Diff{§^). The Lie-group exponential map v — > expj;^(t') is a smooth map 
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of the Lie algebra to the Lie group 123'. Although the derivative of exp^;^ at 
G C°°(S^) is the identity, exp^^ is not locally surjective j23j. This failure, in 
contrast with the case of Hilbert Lie groups ^Hl; is due to the fact that the 
inverse function theorem does not necessarily hold in Frechet spaces jl5j . 

3.2. H'' metrics on Diff{S^). For > and u,v e Vect{S^) = C°°{E,^), 
we define 

k 

(25) {u,v)k= f V(9>)(9>)dx= / Ak{u)vdx, 
where 

(26) A, = 1-A_ + ... + (-!)'=" 



is a continuous linear isomorphism of C°°(S^). Note that is a symmetric 
operator for the inner product 

Ak{u)vdx= I uAk{v)dx. 

Remark. What should be g* for G = Diff{S^) and g = vect{S^) ? If we let 
g* be the space of distributions, A}^ is no longer an isomorphism. This is the 
reason why we restrict g* to the range of A^^ 

Im{Ak) = C°°(§i). 

The pairing between g and g* is then given by the inner product 

(m, u) = / mu dx. 

With these definitions, the coadjoint action of Dif f{§^) on g* = C°°(S^) is 
given by 

^^>= (c^.(^l))^ "^^'^''^- 
Notice that this formula corresponds exactly to the action of the diffeomor- 
phism group Di//(S^) on quadratic differentials of the circle (expressions 
of the form m{x)dx'^). This is the reason why one generally speaks of the 
torsor space of the group Dif f{^^) as the space of quadratic differentials. 


We obtain a smooth right-invariant metric on Dif f{^^) by extending the 
inner product (|25j) to each tangent space T^p Dif f^B^), ip € Diff{§^), by 
right-translations i.e. 

{V, W)^ = {R^-iV, R^-i W)^, V,WG TpDiffiS^). 

The existence of a connection compatible with the metric is ensured (see 
dni) by the existence of a bilinear operator B : C°°(Si) x C°°{S^) C°^{S^) 
such that 

{B{u,v), w) = {u, [v,w]), u,v,we yect(S^) = C7°°(S^). 
For the metric, this operator is given by (see 

(27) Bk{u,v) = -A^^(2v^Akiu)+vAk{u^)y u,vgC'^{S^). 
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3.3. Geodesies. The existence of the connection V'^ enables us to define 
the geodesic flow. A C^-curve ip : I ^ Dif f{^^) such that = 0, where 

denotes the time derivative ipt of (/?, is cahed a geodesic. As we did in 
Section |21 in the case of a left-invariant metric, we let 

u{t) = R^-i ip = ifto Lp-^ 

which is the right angular velocity on the group Diff(S^). Therefore, a 
curve ip G C^(/, Diff{§^)) with p{0) = Id is a geodesic if and only if 

(28) ut = Bk{u,u), t€l. 

Equation (|28j) is the Euler-Arnold equation associated to the right-invariant 
metric (|25]). Here are two examples of problems of type (|28]) on Diff{E>^) 
which arise in mechanics. 

Example. For k = 0, that is for the right-invariant metric, equation ()28() 
becomes the inviscid Burgers equation 

(29) ut + 3Mnx. = 0. 

All solutions of (|29j) but the constant functions have a finite life span and 
(|29)1 is a simplified model for the occurrence of shock waves in gas dynamics 
(see [Ej). 

Example. For k = 1, that is for the right-invariant metric, equation ()28|) 
becomes the Camassa-Holm equation (cf. (21]) 

(30) Ut + uur, 3^ (1 - d^2y'^ + ^ u^2^ = 0. 

Equation (|3U() is a model for the unidirectional propagation of shallow wa- 
ter waves jSjilTj. It has a bi-Hamiltonian structure ,14^ and is completely 
integrable . Some solutions of (|5n|) exist globally in time [3 [H] , whereas 
others develop singularities in finite time [HI El HI HH- The blowup phe- 
nomenon can be interpreted as a simplified model for wave breaking - the 
solution (representing the water's surface) stays bounded while its slope 
becomes unbounded [S]. 

3.4. The momentum. As a consequence of the right-invariance of the met- 
ric by the action of the group on itself, we obtain the conservation of the 
left angular momentum tul along a geodesic p. Since tul = Ad*^-i f^R and 
ruR = Afc(n), we get that 

(31) mk{p,t) = Ak{u) o p ■ p,l, 
satisfies rak{t) = mfc(O) as long as mk[t) is defined. 

3.5. Existence of the geodesies. In a local chart the geodesic equation 
can be expressed as the Cauchy problem 

(32) I = p , 

with p{0) = Id, v{0) = n(0). However, the local existence theorem for 
differential equations with smooth right-hand side, valid for Hilbert spaces 
[THj . does not hold in C°°{S^) (see ^S]) and we cannot conclude at this 
stage. However, in [T], we proved 



LIE GROUPS AND MECHANICS 



15 



Theorem 3.1. Let k > 1. For every uq € C°°(S^), there exists a unique 
geodesic ip G C°°([0, T), for the metric starting at ip{0) = 

Id E Diff{S^) in the direction uq = ipt{0) G Vect{S^). Moreover, the 
solution depends smoothly on the initial data uq £ C^{S^). 

Sketch of proof. The operator in ()32() is specified by 

Pk{^,v) = Qk{v o (^^1) o if, 

where Qk : C°°(S^) ^ C°°(S^) is defined by Qkiw) = Bk{w,w) + ww^. 
Since 

c°°(si) = Pi i/-(si) 

r>n 

for all n > 0, we may consider the problem H32|) on each Hilbert space 
If A; > 1 and n > 3, then is a smooth map from x H'^{S^) 
to H^(S^), where C //"(S^) is the open subset of all functions having 
a strictly positive derivative. The classical Cauchy-Lipschitz theorem in 
Hilbert spaces ^H] yields the existence of a unique solution if nit) G C/" of 
(|32j) for all t € [0, T^) for some maximal T„ > 0. Relation can then be 
used to prove that Tn = T^+i for all n > 3. □ 

Remark. For A: = 0, in problem (|32() . we obtain 

which is not an operator from [/" x H^[S^) into //"(S^) and the proof 
of Theorem 13.11 is no longer valid. However, in that case, the method of 
characteristics can be used to show that even for A; = the geodesies exists 
and are smooth (see JHI)- 

3.6. The exponential map. The previous results enable us to define the 
Riemannian exponential map eyp for the H'^ right-invariant metric (k > 0). 
In fact, there exists 6 > and T > so that for all uq € Diff(§i^) with 
ll^olbfc+i < 5 the geodesic ip{t;uo) is defined on [0,T] and we can define 
e^p{uo) = ip{l;uo) on the open set 

= |no E DiffiS') : holbfe+i < 

of Diff{S^). The map uq i— s- ej:p(no) is smooth and its Frechet derivative 
at zero, DeppQ, is the identity operator. On a Frechet manifold, these facts 
alone do not necessarily ensure that epp is a smooth local diffeomorphism 
jl5j . However, in jJJ, we proved 

Theorem 3.2. The Riemannian exponential map for the right-invariant 
metric on DiffCB^), k > 1, is a smooth local diffeomorphism from a neigh- 
borhood of zero on Vect{§^) to a neighborhood of Id on Diff(S^). 

Sketch of proof. Working in H^^^{S^), we deduce from the inverse function 
theorem in Hilbert spaces that ej:p is a smooth diffeomorphism from an 
open neighborhood Oyt+a of G to an open neighborhood @k+3 

of Id G 

We may choose O^+s such that De^p^^ is a bijection of for 
every uq G Ofc+s- Given n > A; + 3, using (|ST|) and the geodesic equation. 
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we conclude that there is no uq G iJ"(Si) \ iJ"+i(Si), with eyp(iio) G ?7"+^ 
We have proved that for every n > k + 3, 

zTP :0 = Ok+3 n c°°(si) ^ e = e^+g n c°^i§') 

is a bijection. Using similar arguments, (|^T|) and the geodesic equation can 
be used to prove that there is no uq E H"{S^)\H'^+'^{E,^), with De^p^^{v) G 
^71+1 ^gi^ for some lio € O. Hence, for every € O and n > k + 3, the 
bounded linear operator De-^p^^ is a bijection from ff"(S^) to □ 

Remark. For A; = we have that ej:p is not a local diffeomorphism from 
a neig hborhood of € Vect(§^) to a neighborhood of /d € Diff{E^), as 
proved in [lOj. The crucial difference with the case {k > 1) lies in the fact 
that the inverse of the operator A^, defined by H26|) . is not regularizing. 
This feature makes the previous approach inapplicable but the existence of 
geodesies can nevertheless be proved by the method of characteristics. 
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